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Schistosoma mansoni is a parasitic plathyhelminth responsible for intestinal 
schistosomiasis (or bilharzia), a disease affecting 67 million people worldwide and 
causing an important economic burden. The schistosomicides hycanthone, and its 
later proxy oxamniquine, were widely used for treatments in endemic areas during the 
twentieth century. Recently, the mechanism of action, as well as the genetic origin of a 
stably and Mendelian inherited resistance for both drugs was elucidated in two strains. 
However, several observations suggested early on that alternative mechanisms might 
exist, by which resistance could be induced for these two drugs in sensitive lines of 
schistosomes. This induced resistance appeared rapidly, within the first generation, 
but was metastable (not stably inherited). Epigenetic inheritance could explain such a 
phenomenon and we therefore re-analyzed the historical data with our current knowledge 
of epigenetics. In addition, we performed new experiments such as ChlP-seq on 
hycanthone treated worms. We found distinct chromatin structure changes between 
sensitive worms and induced resistant worms from the same strain. No specific pathway 
was discovered, but genes in which chromatin structure modifications were observed are 
mostly associated with transport and catabolism, which makes sense in the context of the 
elimination of the drug. Specific differences were observed in the repetitive compartment 
of the genome. We finally describe what types of experiments are needed to understand 
the complexity of heritability that can be based on genetic and/or epigenetic mechanisms 
for drug resistance in schistosomes. 
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INTRODUCTION 

Schistosoma mansoni, a platyhelminth trematoda parasite, is the 
causing agent of intestinal schistosomiasis (or bilharzia). The par- 
asite has a heteroxenous life cycle with two hosts and three main 
developmental stages. The cycle begins when eggs come in con- 
tact with fresh water and hatch, releasing a free-swimming larva 
known as miracidium. This larva infects the intermediate host 
(a freshwater snail, usually of the Biomphalaria genus). Inside 
the mollusk, it will transform into primary and then secondary 
sporocysts and undergo asexual multiplication. A second type 
of free-swimming larvae, cercariae, is released and infects the 
definitive mammalian host (human or rodent) in which they 
will reach their adult form (worm) and perform sexual repro- 
duction. Sexually mature schistosomes couples are located in 
mesenteric veins of the mammalian host. Females lay eggs with 
a characteristic lateral spine, which allow eggs to go through the 
wall of the intestine and be excreted in the feces. However, a 



large proportion of eggs is caught in the blood circulation and 
end up being trapped in the liver, leading to hepatic fibrosis. 
Bilharzia is an endemic disease in many areas in Africa and South 
America, infecting 67 million people and causing an impor- 
tant socio-economical burden (King, 2010). Since the beginning 
of the twentieth century, several schistosomicides were devel- 
oped (Galdino da Rocha Pitta et al., 2013). Many of them were 
abandoned due to low treatment efficiency, as well as toxic and 
sometimes mutagenic side effects to the human host (Galdino 
da Rocha Pitta et al, 2013). Another major concern regarding 
schistomicides is the rapid appearance of resistances, both in lab- 
oratory strains and natural populations of the parasite (Bruce 
et al, 1987; Coles et al, 1987; Drescher et al, 1993; WHO, 1998). 
Please note that to avoid any language confusion, we use here 
the definitions of drug resistance and tolerance in schistosomes as 
given by Coles (2002), with the former being "when any isolate has 
a significantly lower cure rate than the most sensitive isolate" and the 
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later when the whole species is non-responsive to the drug. These 
factors left a single drug, praziquantel, as the current therapeutic 
alternative. As there is very little advance in the development of 
novel drugs against schistosomiasis, it becomes crucial to deepen 
our understanding on which molecular and evolutionary mech- 
anisms can lead to schistosomicides resistance. Such knowledge 
could allow to either avoid the rapid appearance of resistances 
previously seen in S. mansoni, or to find ways to circumvent it. 

Various drug resistances in schistosomes were studied through 
the past decades, but one of the best documented, but also most 
intriguing one, is against hycanthone. This molecule is not used 
anymore in treatments, and has been commercially discontinued, 
but was subjected to many research works from different labo- 
ratories. First evidences of resistant worms to this molecule were 
published in the early 1970's (Rogers and Bueding, 1971). It was 
later demonstrated several times that schistosomes were in gen- 
eral cross-resistant or cross-susceptible to both hycanthone and 
another schistosomicide, oxamniquine (Jansma et al., 1977; Cioli 
and Pica-Mattoccia, 1984; Dias and Olivier, 1985; Bruce et al., 
1987). Oxamniquine was often used as a proxy for hycanthone, 
as the later is not commercially produced anymore and was toxic 
to the host to a certain degree. Two laboratories obtained stable 
resistant lines of schistosomes by injecting infected mice with a 
single curative dose (i.e., killing >90% of the worms in suscepti- 
ble lines, 80 mg/kg) repeated over 3-5 successive generations of 
the parasite (Cioli and Pica-Mattoccia, 1984; Dias and Olivier, 
1985). In these lines, resistance was heritable and stable even in 
the absence of drug pressure over 9-30 generations (Dias and 
Olivier, 1985; Cioli et al, 1989; Drescher et al, 1993). Classic 
genetic crosses between sensitive and resistant strains, as well 
as genetic complementation demonstrated that a single recessive 
autosomal locus codes for the hycanthone/oxamniquine resis- 
tance in two strains, MAP (named from the initials of the patient 
from which it was isolated) and Baltimore Rome Resistant (BRR) 
(Cioli and Pica-Mattoccia, 1984; Cioli et al, 1989, 1992; Pica- 
Mattoccia et al., 1992b, 1993). Similar results were also obtained 
in vitro by Coles and Bruce (1987). They cultured schistosomulas 
(an immature developmental stage occurring once cercariae have 
entered the mammalian host) under 240 p.p.m. of oxamniquine, 
a lethal concentration for 99.8% of the population. Survivors 
were injected in mice, and their progeny was not responsive to 
very high concentrations (500 mg/kg) of the drug. Experiments 
with tritiated hycanthone and oxamniquine demonstrated that 
the molecules were absorbed by all worms, but did not persist in 
resistant worms, while they were found to be covalently bound to 
DNA, RNA, and proteins in the sensitive strains (Pica-Mattoccia 
et al., 1988, 1989). Biochemical characterization showed that 
an enzymatic activation, probably by a sulfotransferase (Pica- 
Mattoccia et al., 2006), was needed by both oxamniquine and 
hycanthone to be lethal to the worm. Prior to its formal character- 
ization, Pica-Mattoccia et al. (1992a) were able to biochemically 
isolate a fraction containing this enzyme from ground sensi- 
tive worms. They showed that resistant worms were killed when 
exposed in vitro to either hycanthone or oxamniquine and to the 
sulfotransferase-containing fraction from sensitive worms, while 
unaffected when exposed to one of the drugs alone. In the same 
article, the authors established that adult male worms, that are 



more sensitive to the drugs than females, have a higher enzymatic 
activity than females and immature stages (mostly unaffected by 
the drugs at the classical therapeutic concentration, and having 
a very low enzymatic activity). A difference in the enzyme speci- 
ficity might also explain why hycanthone, but not oxamniquine, is 
efficient against Schistosoma haematobium (Pica-Mattoccia et al., 
1997). The sequencing of the genome of S. mansoni (Berriman 
et al., 2009) provided new opportunities of investigation, and 
Valentim et al. (2013) recently used linkage mapping methods to 
identify the gene responsible for resistance in the HR and MAP 
strains. Candidate genes were screened using recombinant pro- 
teins, RNAi, and protein crystallography. The authors showed 
that it was indeed a sulfotransferase that was responsible for the 
resistance. They identified its coding gene (Smp_089320) and 
the non-synonymous nucleotide substitutions causing the loss 
of function of the enzyme in resistant lines MAP and HR, thus 
clearly establishing the mechanism of resistance in the studied 
lines (see Figure SI in Valentim et al., 2013 for the mode of 
activation of oxamniquine). 

Nevertheless, there are several observations that suggested 
early that alternative mechanisms might exist. For instance, 
Jansma et al. (1977) described three different ways to induce resis- 
tance by collecting the progeny of worms after a single injection of 
60 mg/kg of hycanthone to the rodent host 54-70 days post infec- 
tion (Type I), or 27-29 days post infection when worms are still at 
an immature stage (Type II), and from hosts infected by a single 
sex and then re-infected 2-58 weeks later by cercariae of the oppo- 
site sex (Type III). In all three cases, parents were all susceptible 
to the drug while a high percentage of the progeny survived the 
administration of a curative dose of hycanthone. Resistance was 
shown to be metastable to various degrees (depending on the type 
of induction used) through 10-21 generations. By metastable, 
we mean that from one generation to another, the percentage 
of worms surviving hycanthone treatment was highly fluctuant 
and did not follow a classical Mendelian inheritance pattern. For 
example, it was observed at several occasions by Jansma et al. 
(1977) that while a parental generation was completely resistant 
(100% survivors), a high proportion of the direct progeny (some- 
times up to 75%) was found sensitive. Sometimes, the amount 
of resistant parasites would go up again a couple of generations 
later. These results are in complete contrast of the work men- 
tioned above where classical recessive monogenic pattern for the 
inheritance of the resistance was observed. The work of Jansma 
was partly reproduced by Brindley et al. (1989). They were able 
to achieve Type II induction with lines of schistosomes deriving 
from Jansma's isolate. 

As these observations do not fit the monogenic resistance 
model recently described, we reasoned that environmentally trig- 
gered drug resistance with metastable inheritance is a legitimate 
candidate for an epigenetically based phenotype. We use here a 
more narrow definition of epigenetics, meaning heritable changes 
in gene expression, without changes in the DNA sequence (Wu 
and Morris, 2001). They can be physically based on chemical 
modification of the DNA (e.g., methylation of cytosines), changes 
of proteins interacting with DNA (e.g., histone modifications), as 
well as nuclear localization of chromosomes, and short untrans- 
lated RNA playing a role in post-transcriptional silencing of 
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genes (Aguilera et al., 2010). There are numerous examples in 
the literature on how environmental factors, such as nutrition, 
pollution, temperature, and others, can have an influence on epi- 
genetic control (Feil and Fraga, 2012). A drug treatment could, 
for example, induce an epimutation, silencing the gene coding for 
the enzyme responsible of the bio-activation. This epimutation 
would be more or less stable and/or inherited to progeny given 
the parasite strain (as the phenotype was previously described 
to be either stable or metastable through generations depending 
on investigators and strains). This type of transient transgenera- 
tional modification was initially also called "Dauermodifikation" 
(Jollos, 1921). 

From this hypothesis, we decided to revisit the historical data 
and to complement it with our own study. DNA methylation has 
been found to regulate oviposition (Geyer et al., 2011). However, 
it is globally at a very low level through the genome (Raddatz 
et al., 2013), and therefore we decided to focus on another epi- 
genetic information carrier, histone modifications. We decided to 
use a Brazilian (SmBre) and Guadeloupian (SmGH2) strain. Both 
strains were previously pool-sequenced (Clement etal., 2013) and 
we found that they were both of low genetic diversity (Tajima's pi 
of 2.10" 4 for SmGH2 and 1.8. 10" 4 for SmBRE) and do not have 
any coding polymorphism in the gene (Smp_089320) uncovered 
by Valentim et al. (2013). Even if there is an almost uniform 
genetic background within the strains, we could not, at this stage 
completely rule out that an induction phenotype could come 
from the selection of a novel resistance allele in some schistosomes 
and not an epigenetic change. The debate whether heritable char- 
acters are induced by the environment or whether they occur 
randomly and are subsequently selected is not new. In 1951, 
in order to differentiate between "spontaneous (genetic) muta- 
tion and natural selection" and "directed mutation" or induction, 
the Lederbergs had designed their now famous experiment using 
replica plating of bacterial colonies from non-selective support 
to Petri dishes with selective media (Lederberg and Lederberg, 
1952). Their work gave results in favor of the spontaneous muta- 
tion, i.e., a genetic variability (diversity) based mechanism. We 
reasoned that it would be in principle possible to reproduce 
the classical Lederberg experiment, but this time using clonal 
populations of S. mansoni instead of bacterial clonal colonies. 

Our results show that one of the strains, SmBre, is resistant 
to hycanthone in the absence of mutations in Smp_089320. After 
exposure to the drug, the chromatin structure is altered around 
genes belonging to gene ontology (GO) groups involved in cat- 
alytic activity and response to stress. This finding suggests that 
epigenetic mechanisms could be responsible for adjustments of 
the parasites metabolism leading to temporary resistance, thus 
reconciliating the seeming contradictions between stable and 
metastable types of the phenotype. 

MATERIALS AND METHODS 
MAINTENANCE OF THE PARASITES 

The Schistosoma mansoni strains SmBre and SmGH2 were origi- 
nally sampled in Brazil and Guadeloupe, respectively. Both strains 
were maintained in their sympatric intermediate host strain 
BgBre and BgGua of the mollusk Biomphalaria glabrata and 
Mus musculus as definitive vertebrate host. SmBre was originally 



sampled in the 1960s in Recife, and SmGH2 was isolated on Nov. 
10th, 1983 from a patient in Ste. Rose (Guadeloupe, French West 
Indies) (Pers. Commun. A. Theron). 

PRELIMINARY DRUG SCREENING OF THE STRAINS 

Both SmBre and SmGua were screened for resistant individuals 
by first infecting five mice (for each strain) with 70 mixed-sex 
cercariae. Mollusks from which cercariae were collected were ini- 
tially infected with 20 miracidia each, and hence cercariae shed 
are of multiple genotypes. Mice were injected 62 days later with 
a dose of 60mg/kg of hycanthone (similar to the one used by 
Jansma et al, 1977). Eight weeks later, the mice were sacrificed 
by a lethal intraperitoneal injection of sodium pentobarbital, and 
surviving adult worms were recovered by retrograde perfusions of 
the hepatic portal system with citrate (7.5%) saline (8.5%) solu- 
tion administrated through the left ventricle (Duvall and DeWitt, 
1967). Worms trapped in the liver or mesenteric system were 
collected after excising these organs. 

LEDERBERG-LIKE EXPERIMENT 

Biomphalaria glabrata BgBre individuals were infected with single 
miracidium from SmBre. The experiment was repeated with sev- 
eral different B. glabrata BgBre host individuals, each one infected 
with a different S. mansoni clone (each miracidium came from 
different parents). We genotyped with sex markers cercariae shed 
by each B. glabrata following Beltran et al. (2008), and we selected 
five mollusks infected with male parasites. We chose to work only 
on male schistosomes as it has been shown that they are more sen- 
sitive to hycanthone (Pica-Mattoccia et al, 1992a) than females, 
and to avoid sex-based bias in the epigenome analysis. From each 
snail, we infected 10 mice with 70 cercariae each. Sixty-two days 
later, half of the mice were injected with a dose of 60 mg/kg of 
hycanthone, and the other half with Ringer's solution as a nega- 
tive control. Fifty-seven days after injection, we perfused the mice 
following the method of Duvall and DeWitt (1967) to collect con- 
trol and surviving, resistant worms. A graphical abstract of the 
methodology is on Figure 1. 

CHROMATIN IMMUN0PRECIPITATI0N AND SEQUENCING 

Native ChIP was performed following the protocol developed 
for Schistosoma mansoni by Cosseau and Grunau (2011) (also 
available online at http://methdb.univ-perp.fr/cgrunau/methods/ 
native_chip_sm.html) on a pool of six resistant and a pool of 
10 control adult worms. Immunoprecipitation was performed 
using the following antibodies: H3K4me3 (Millipore, cat# 04-745 
lot* NG1680351, 4 ul per reaction), H3K9Ac (Millipore, cat# 07- 
352 lot* DAM16924924, 8 ul per reaction), H3K9me3 (Abeam, 
cat* ab8898 lot* 733951, 4ul per reaction), and H3K27me3 
(Diagenode, cat* pAb-069-050 lot* A29900242, 8 ul per reac- 
tion). ChIP products were sequenced as single-end 50 bp reads 
on an Illumina HiSeq 2000 at the Montpellier GenomiX facil- 
ity (http://mgx.cnrs.fr). We performed PCR amplification of the 
two exons of Smp_089320 on control and resistant ChIP prod- 
ucts to test if a possible somatic mutation had occurred in that 
gene. We used Promega GoTaq Flex (cat* M8291) at 3mM 
MgCl2 and with the following PCR conditions: initial denatura- 
tion 2min at 95° C followed by 35 cycles of denaturation 95° C 
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Mono-miracidial infection 



Identification of sex by PCR, 
Mollusks with male 
sporocysts are kept 

Infection with 70 cercariae 



10 mice 



Injection with 60 mg/kg of hycanthone 
or Ringer's solution (control) after 62 days 



Perfusion to recover surviving 
worms 57 days after injection 



FIGURE 1 | Graphical abstract of the methodology used for the 
Lederberg-like experiment. 



for 30 s, annealing 62° C for 30 s and elongation 72° C for 1 min. 
A final extension of 72° C for 5 min completed the reaction. We 
used two primer pairs to amplify exon 1 (Fwdl tccacctctctctcact- 
caatg, Rev 1 ccacatgttggtaaatccgta, Fwd2 gctggtctaccgagaactgg and 
Rev2 tcccgtagaaaaccaactcg) and one pair for exon 2 (Fwd agtc- 
cattcattcaatgtttcaa and Rev caatccacaaatccccattc). PCR products 
were gel purified using Promega Wizard SV Gel and PCR Clean- 
up System (cat# A9282) and sequenced at the Genoscreen facility 
(http://www.genoscreen.fr/). 

QUALITY CONTROL AND ALIGNEMENT 

All data treatment was carried out under a local galaxy instance 
(Goecks et al., 2010) (http://galaxy.univ-perp.fr). We used the 
FASTX-toolkit (http://hannonlab.cshl.edu/fastx_toolkit/) for ver- 
ification of the reads quality. Read quality was judged sufficiently 
good (the majority of reads showed a fastqsanger quality score 
>24 for all positions). For the analysis of unique sequences, reads 
were aligned on the Schistosoma mansoni genome assembly v5.0 
(Protasio et al., 2012) with Bowtie2 (Langmead and Salzberg, 
2012) using parameters — sensitive -k 2. Reads with unique match 
were filtered with samtools (Li et al., 2009) (samtools view -hS - 
q 255). For peak calling, an equal number of randomly selected 
reads was chosen between the treated and control conditions 
(H3K4me3 = 2.9 million, H3K9Ac = 23 million, H3K9me3 = 
21 million, H3K27me3 = 17.9 million). Peak identification was 
performed using PeakRanger vl.16 (Feng et al., 2011) with P- 
value cut off 0.0001, FDR cut off 0.01, Read extension length 
200, Smoothing bandwidth 99 and Delta 1. Wiggle files generated 
by PeakRanger where uploaded in our local instance of GBrowse 



(Stein et al., 2002) (http://genome.univ-perp.fr) for visualization. 
For the analysis of repetitive sequences, reads were aligned with 
Bowtie2 evoking parameters — sensitive -k 5. The 3,145 repetitive 
consensus sequences used as reference in this study were obtained 
from a previous analysis done with the same strain SmBre 
(Lepesant et al, 2012). From the resulting SAM file, we choose 
the position with the highest score for each read. Normalization 
of SAM files was done by randomly removing a given portion of 
the counts to achieve the same effective library size (H3K4me3 
= 0.7 million, H3K9Ac = 9.8 million, H3K9me3 = 10.8 million, 
H3K27me3 = 13.9 million). Elements that showed differences in 
histone modification levels where annotated using the ABBlast 
search engine within the RepeatMasker program (Smit et al., 
2010). The search was performed for both nucleotide and pro- 
tein databases. Sequences without match where compared against 
the nucleotide database of NCBI using the blastn search engine 
and against the non-redundant protein database using blastx 
(Altschul et al, 1990). 

COMPARATIVE ANALYSIS 

To find differences in the histone mark pattern between con- 
trol and treated conditions, we used the Pyicoenrich function 
(with -binsize 0.3) from Pyicos v2.0.6 (Althammer et al., 2011). 
Visual inspection was performed on all the regions detected by 
the software with a z-score above 7. The GO enrichment anal- 
ysis was performed using BLAST2GO (Conesa et al., 2005). For 
the analysis of significant differences in histone isoform enrich- 
ment in repeats, normalized read counts were compared using the 
DESeq package of the R software (Anders and Huber, 2010). This 
package uses the negative binomial distribution to calculate the P- 
value for each element between two compared samples and find 
elements that present significant differences in their enrichment 
levels. 

GO ENRICHMENT ANALYSIS 

The 60 genes for which chromatin structure differences had been 
detected were annotated using Blast2GO (Conesa et al., 2005). 
In order to check if GO terms would be overrepresented in this 
group of genes, 60 transcripts were randomly sampled out of 
41,669 transcripts generated by a combination of RNA-Seq data 
and Sanger annotation. RNA-Seq data were generated from sev- 
eral biological replicates and these data will become available 
when the analysis is fully finished. The random samples were 
annotated as above and GO enrichment analysis was performed 
using Fisher's exact test implemented in Blast2GO and a p-value 
of 0.2 as cutoff. Random sampling was repeated 3 times. 

RESULTS 

A STRAIN WITH SENSITIVE GENOTYPE DISPLAYS HYCANTHONE 
RESISTANCE 

We compared the number of worms recovered in mice injected 
with hycanthone to those recovered from control mice. We found 
25 (36%) surviving worms, perfused from 5 mice, in the Brazilian 
strain SmBre, despite the fact that it has the sensitive genotype 
for Smp_089320. The Guadeloupian strain SmGH2, however, is 
completely sensitive to hycanthone, i.e., no survivors were found 
in any of the five infected and treated mice. 
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CLONAL SCHISTOSOMES RESPOND DIFFERENTLY TO HYCANTHONE 
TREATMENT 

Following the discovery of hycanthone resistant individuals 
within the SmBre population, we adapted the Lederbergs' exper- 
iment to our schistosome model. To do so, we infected five 
Biomphalaria glabrata with a single, genetically distinct, SmBre 
miracidia. Cercariae shed by each mollusk are genetically iden- 
tical. We then infected 10 mice with 70 cercariae. Half of the 
mice were then treated with an injection of hycanthone similar 
to the one described by Jansma for Type I induction (Jansma 
et al., 1977), and the other half was injected with Ringer's solu- 
tion, as negative controls. As we excluded the possibility of a 
somatic mutation in Smp_089320, similar to the classical replica 
plating experiment, two possible outcomes can be imagined: (1) 
if a random mutation had occurred in one of the clones that 
conferred resistance to hycanthone (pre-existence), then only cer- 
cariae from this very snail would produce resistant flukes because 
they are genetically identical (Figure 2A) (2) if the phenotype 
is induced through the treatment with a certain probability, we 
would observe resistance randomly associated with different mice 
hosts and without correlation with the particular snail host (sim- 
ilar to Lederbergs' bacterial colony, Figure 2B). We later perfused 
the mice to look for control and surviving, resistant worms. The 
results of our experiment are shown in Figure 3. In summary, 
while the vast majority of mice were exempt of parasites after 
hycanthone injection, we found surviving worms in six different 
rodents, belonging to three different S. mansoni genotypes (B- 
D). We did not find any worm originating from genotype A and E 
after the treatment. In contrast, under the control condition, we 
recovered an average of 30 ± 5 worms (confidence interval at a 
5%) per mice for the five genotypes. Sequencing of Smp_089320 
in resistant worms did not show any of the three mutations 
identified by Valentim et al. (2013) or any new polymorphism. 

CHROMATIN STRUCTURE CHANGES ASSOCIATED WITH THE 
RESISTANT PHENOTYPE 

We performed ChlP-Seq on the pool of the six hycanthone 
resistant worms, as well as on control worms. We targeted 



four modifications of histone H3: tri-methylation of lysine 4 
[H3K4me3], tri-acetylation of lysine 9 [H3K9ac], tri-methylation 
of lysine 9 [H3K9me3], and tri-methylation of lysine 27 
[H3K27me3]. The first two are associated with euchromatin 
(relaxed DNA structure allowing transcription) and the later two 
are associated with heterochromatin (condensed DNA structure 
impeding transcription). Around Smp_089320 none of the marks 
showed enrichment. We then searched for variations in the chro- 
matin structure at the genome level, and a bioinformatic analysis, 
followed by visual inspection, was able to detect differences in 
the histone modification profiles between hycanthone resistant 
and control worms. We found 47 differences for H3K4me3, 43 
of which being located at the transcription start site of genes. 
This corresponds roughly to 0.6% of the total peaks detected 
for this histone mark. For H3K9me3, 12 differences were spot- 
ted (0.012% of total peaks), mostly in introns and within 2kb 
upstream and downstream of genes. We identified five differences 
for each of the other two marks (0.007% and 0.005% of total 
peaks for H3K9Ac and H3K27me3, respectively), the majority 
being outside of genes. A list of the regions and associated genes 
in which we found differences is in Supplementary Table 1. 

CHROMATIN STRUCTURE MODIFICATIONS ARE PREFERENTIALLY 
ASSOCIATED WITH GENES RELATED TO CATALYTIC ACTIVITY 

None of the genes in which differences were uncovered code for a 
sulfotransferase. Using a GO enrichment analysis, we did not find 
any clear pathway linking together the genes bearing chromatin 
structure changes. However, we found that the terms "transferase 
activity" (13 occurrences), "kinase activity" (8 occurrences) and 
"catalytic activity" (26 occurrences) and "response to stress" (5 
occurrences) were overrepresented in the panel of genes with 
epimutations. 

CHROMATIN STRUCTURE MODIFICATIONS IN REPETITIVE SEQUENCES 

We identified, in 14 repetitive sequences, statistically signifi- 
cant changes in chromatin structure. This is very few (0.44% 
of repeats) showing that globally the chromatin remains stable. 
Interestingly, changes in H3K4me3 and H3K9ac always decrease 
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FIGURE 2 | Schematic representation of the experimental strategy for 
study of the origin of hycanthone resistance. Hypothesis (A) rare 
mutations (star) in miracidia (top) confer resistance to hycanthone. All 
offspring will be resistant (or sensitive, not shown). Hypothesis (B) 
epimutations or mutations occur during clonal amplification of cercaria in the 



snail host. Only few adult worms will be resistant. In the case of epimutation, 
in subsequent generations reversal of the phenotype would occur (and was 
described in the past) and a chromatin structure changes can be detected. In 
case of mutation, no reversal occurs and mutations should be detected in the 
locus of interest. 
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FIGURE 3 1 Results of the experiment on induction of hycanthone 
resistance using mono-miracidial infections. Five mollusks were 
infected using a single miracidium (genotypes A-E). Cercariae shed 
from a single mollusk are clonal and were used to infect 10 mice. 
Half of the mice were injected with a control Ringer's solution and 



the other half with a therapeutic dose of hycanthone. Mice were 
later perfused and worms were counted. In the hycanthone 
condition, six worms of three genotypes (B-D) were recovered from 
six mice, while an average of 30 worms per mouse were found in 
the control condition. 



in the treated worms while H3K27me3 increases, indicating a 
tendency toward heterochromatization. No changes occur in 
H3K9me3. H3K4me/H3K9ac modifications are observed in dis- 
persed repeats that belong to the LINE/LTR classes. In contrast, 
H3K27me3 becomes enriched in simple repeats that are orga- 
nized in blocks (Supplementary Table 2) sometimes covering the 
entire length of large contigs for which chromosomal location is 
unknown. 

DISCUSSION 

When screening the literature on hycanthone resistance, we were 
surprised to see such divergent results between the different 
research groups. On one hand, there is clearly established mono- 
genic, stably inherited, resistance, and on the other, there is an 
inducible resistance, appearing from sensitive parents, and not 
inherited in a Mendelian fashion. We took a deeper look at the 
methodology employed by the all the research groups, and their 
interpretation of data, with the hope of finding variations in 
the procedures, or in the analysis methods, that could, at least 
partly, explain the discrepancies between the observations and 
conclusions. 

We first tried to track back the origins of the strains that were 
employed in the various studies on hycanthone/oxamniquine 
resistance (Table 1). Our conclusion is that very few strains were 
used in these studies, and most of them had been maintained in 
laboratory for a very long time, probably reducing the genetic 
diversity of the strains. We identified six resistant strains in the 
literature, and four of them derived from theM strain used in 
Jansma et al. (1977) for the first induction experiments. We also 
found four main sensitive strains. At least four other strains 
(either resistant or sensitive) were punctually used and will not 
be discussed here (Bruce et al., 1987; Drescher et al., 1993; 
Pica-Mattoccia et al, 1993). 



In his Type I induction, Jansma (Jansma et al., 1977) injected a 
dose of 60 mg/kg of hycanthone 54-70 days post-cercarial infec- 
tion. He noted a marked reduction in the number of worms, 
but rarely a complete elimination. Surviving worms were unable 
to produce viable eggs for 4-8 months. After that period, eggs 
from which miracidia were hatched were collected to produce a 
Fl. This first generation was treated with a therapeutic dose of 
hycanthone (80 mg/kg) and 94% of the worms were found to 
be resistant. Twenty-one successive generations were produced 
from this Fl, and on several occasions, additional groups of mice 
were infected from cercariae batches shed on different days, by 
the same infected mollusks, leading to parallel generations known 
as branches. Resistance tends to decrease with generations in the 
main branch, and is highly variable between branches (some are 
totally resistant for several generations, some become completely 
sensitive, and others are fluctuant). Although it is not explicitly 
stated in the paper, it seems that progeny for generations F2-F2 1 
was collected before the therapeutic drug injection used to esti- 
mate the percentage of resistant parent worms. Hence, there was 
no drug selective pressure from one generation to another, and 
that could partly explain the metastability observed by Jansma 
but not by other investigators, as suggested by Cioli and Pica- 
Mattoccia (1984). It is conceivable that Type I resistance was 
achieved not by induction, but through the selection worms with 
a resistant genotype, and as the induction dose is slightly below 
the therapeutic one (normally killing >90% of the parasites), 
there is possibility that some worms with a normally "sensitive" 
phenotype, escaped the drug. Cioli et al. (1992) observed that 
their resistant strain derived from Jansma's had a lower global 
infectivity (to mollusk and mammalian host) and a lower fertil- 
ity (lower amount of eggs per females) than their sensitive strain 
(although this could be strain related, as the sensitive strain from 
Cioli is not from the same origin as Jansma's). In the hypothesis 
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Table 1 | Origin and resistance phenotypes of the most commonly used strains in the literature. 



Strain 


Origin 


Resistance 


Derived from? 


Comment 


References 


M 


Puerto Rico 


Oxa/Hyc 


- 


Used for the three types 
of induction experiments 


Jansma et al., 1977 


BRR 


Puerto Rico 


Oxa/Hyc 


M 


Underwent drug 
selective pressure over 3 
generations 


Cioli et al., 1992; Pica-Mattoccia et al., 1992a, b, 
1993, 2006 


JHU 


Puerto Rico 


Oxa/Hyc 


M 


Successful Type II 
induction 


Brindley and Sher, 1987; Brindley et a!., 1989, 
1991 


H-30 


Puerto Rico 


Oxa/Hyc 


M 




Souza et al., 1988 


HR 


Puerto Rico 


Oxa/Hyc 


M 




Valentim et al., 2013 


MAP 


Brazil, Minas Gerai 


Oxa Partial to Hyc 






Dias and Olivier, 1985; Bruce et al., 1987; 
Souza et al., 1988; Pica-Mattoccia et al., 
1992a,b, 1993; Drescher et al., 1993; Valentim 
etal., 2013 


SEN 


Puerto Rico 


None 






Cioli, 1976 


BH 


Brazil, Minas Gerai 


None 




Isolated from an 
untreated patient in 1967 


Dias and Olivier, 1986; Bruce et al., 1987 


NMRI 


Puerto Rico 


None 




Isolated in the 1940' 


Fletcher et al., 1981 ; Brindley and Sher, 1987; 
Brindley etal., 1989 


LE 


Brazil 


None 






Valentim et al., 2013 



Hyc, Hycanthone resistant; Oxa, Oxamniquine resistant. 



that some "sensitive" worms survived the induction step and were 
mixed with selected resistant worms, this could explain why the 
resistance is does not follow a Mendelian inheritance. This can- 
not be the case in Type II, as immature worms have little or no 
sensitivity to hycanthone (induction was also successful with sub- 
curative doses ranging from 3 to 60 mg/kg) and in Type III where 
the parents were not exposed to the drug before collecting the 
progeny. These observations lead us to think that the resistant 
phenotype in Type II and III could indeed be induced and trans- 
mitted to the progeny to a certain extent. Cioli and Pica-Mattoccia 
(1984), Cioli et al. (1989) and Dias and Olivier (1986) stated that 
they were unable to reproduce lansma's results in inducing Type 
II or III resistance. Three other antihelminthics, oxamniquine, 
praziquantel, and oltipraz were also unsuccessfully tested with 
the hope of inducing Type II resistance and selecting a new line 
of resistant schistosomes (Dias and Olivier, 1986). However, we 
have a slightly different interpretation of the results from Dias 
and Oliver. In their induction protocol, mice were treated with 
subcurative dose of one of the three drugs, 26-30 days after cer- 
carial infection (when schistosomula have a low sensitivity to both 
hycanthone and oxamniquine), and then divided into two groups. 
The first group was treated with a curative dose of drug, while 
the second was left as a control, and the eggs from that control 
group were used to restart a cycle of the parasite. Mice infected by 
the Fl cercariae from the control group were also divided into a 
trial group treated with a curative dose of the drug and a con- 
trol group (which would serve to generate a F2). No resistant 
worms were found in the treated Fl and F2, but 8-12% where 
found in the parental generation (in which the induction step was 
performed). Admittedly, there was no inheritance of a resistance 
phenotype to Fl and F2, but the fact that survivors were found 
only in the parental generation which underwent the induction 



step let us believe that it was not as a failure at inducing resistance 
as described by the authors. In fact, only when induction was 
performed, some resistant worms were found. Type II was also 
successfully induced by the Brindley group (Brindley et al., 1989) 
in several lines of the parasite deriving from Jansma's isolate. 
However, while four strains tested in this article were described as 
sensitive by the authors, it is worth to mention that in some of the 
trials, up to 26% of the worms survived the hycanthone treatment 
(with or without preliminary Type II induction). Brindley also 
observed spontaneous occurrence of resistance in other strains 
with no induction, although he admitted that "whether these fluc- 
tuations reflect true changes in the parasite genotype rather than 
experimental variation was not systematically investigated in this 
study" (Brindley et al, 1989). This lead us to realize that for the 
majority of research groups, resistance was considered as achieved 
when at least 90% of the worm population survived a therapeu- 
tic dose. In several studies (Dias and Olivier, 1986; Brindley and 
Sher, 1987; Brindley et al, 1989; Drescher et al., 1993), various 
lines of schistosomes did not reach this resistance threshold (with 
or without induction and/or drug pressure over a couple of gener- 
ations). They were considered as sensitive although in some cases 
there was a non-negligible percentage (3-26%) of worms sur- 
viving a therapeutic dose of either hycanthone or oxamniquine, 
strengthening our hypothesis that a mechanism following a non- 
Mendelian inheritance pattern is involved. In this context, we 
believe it is more appropriate to talk about the penetrance of 
the phenotype rather than talking strictly about sensitivity or 
resistance. 

Also, Cioli and Pica-Mattoccia specifically mentioned that 
their approach was different from the Jansma/Bueding/Brindley 
group. Three generations of selection were necessary to obtain 
a phenotype that was stably inherited by the progeny (Cioli and 
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Pica-Mattoccia, 1984). The same observation was made by Dias 
and Olivier (1985), as in their experiments, resistance was sta- 
ble only after the 5th generation (experiments were continued 
until 14th generation). Unfortunately, the data for the first 3-5 
generations are not available but one can assume that incom- 
plete penetrance of the phenotype was observed. These large 
differences in penetrance of resistance phenotypes (i.e., large dose 
dependency) in different isolates were also found by others (Coles 
et al., 1986; Kinoti, 1987) and also described in Pica-Mattoccia 
et al. (1993). Kinoti noted important variations in the dose of 
oxamniquine needed to eliminate 99% of the worms (ED99), 
in both laboratory and clinical trials. For example, East African 
isolates needed a dose of oxamniquine 200-250% higher than 
Puerto-Rican isolates to reach ED99 (Kinoti, 1987). In all publica- 
tions, at least 3 generations of selection were necessary to produce 
stable resistance, but even then, phenotypic differences between 
the different resistance strains and most notably an exception to 
the Mendelian segregation pattern were observed. In the MAP 
strain for instance there was ( 1 ) a difference between heritability 
of hycanthone and oxamniquine resistance and (2) segregation 
did not follow the recessive pattern. Segregation could not be fol- 
lowed up because the females produced no eggs, i.e., no genetic 
analysis was possible ("Thus, we are left with no obvious explana- 
tion for this unexpected partial exception to the overall recessivity 
of the MAP resistance!') (Pica-Mattoccia et al., 1993). The rea- 
son for the incomplete penetrance in the first 3-5 generations 
cannot be determined with confidence from the published data 
alone. It might be that the initial strains showed a high degree of 
genetic diversity with only a few individuals carrying the muta- 
tion in the sulfotransferase gene. Molecular biology methods in 
S. mansoni research were not sufficiently advanced in that time 
to measure genetic diversity. From the few available data, how- 
ever, it seems that the founder populations of the resistant strains 
were very small. Years of maintenance cycle in laboratories, as well 
as the possibility of a founder effect every time that the strains 
were transferred or shared with other laboratories let us think 
that the genetic diversity was probably low in the different strains. 
We also noted that most of the strains were from the same geo- 
graphical areas (either Puerto Rico or the region of Minas Gerais 
in Brazil) and represent only a small share of geographic diver- 
sity of S. mansoni. Another reason for variability in penetrance 
could therefore be high phenotypic plasticity and/or variability in 
plasticity. lansma's results were in favor of this latter mechanism, 
and he additionally described transgenerational plastic effects. It 
has also to be noted that while the work of Brindley and lansma 
was performed using hycanthone, other research groups used 
oxamniquine. Although schistosomes where always cross resis- 
tant or cross sensitive to both molecules, hycanthone is known 
to have a mutagenic effect, and it could have played a role in 
the phenomenon of induction, which was never observed with 
oxamniquine. 

In our own experiments, we worked with a strain bearing 
the sensitive genotype for the sulfotransferase gene previously 
identified as the source of hycanthone/oxamniquine resistance, 
but we still found in drug screening that in one strain (SmBre) 
about a third of the adult worms survived the treatment while 
another strain did not show survivors. Resistance of SmBre 



could be due to a mutation in a different locus and/or due to 
epimutations but this would still not explain Jansma's induc- 
tion phenomenon. Therefore, we used a Lederbergs experiment 
to discriminate if the resistance was pre-existing or induced by 
the experimental conditions. Similar to the classical replica plat- 
ing experiment, two possible outcomes could be imagined: (1) 
if the phenotype is induced through the treatment with a cer- 
tain probability, we would observe resistance randomly associated 
with different mice hosts and without correlation with the par- 
ticular snail host (similar to Lederbergs' bacterial colony) (2) 
if a random mutation had occurred in one of the clones that 
conferred resistance to hycanthone (pre-existence), then only cer- 
cariae from this very snail would produce resistant flukes because 
they are genetically identical. A caveat with the Lederberg exper- 
iment in organisms with epigenetic inheritance systems is the 
potentially heritable transcriptional and/or posttranscriptional 
level of diversity. It is possible that epimutations occur during 
clonal amplification of the cercariae (Figure 2B). Also, as long 
as somatic mutations cannot be excluded, mutation would do 
the same job. Both would mimic the outcome of the "induc- 
tion" alternative. We found that not all the worms issued from 
the same miracidia (i.e., being genetic clones of each other), 
responded to the treatment the same way. While the vast major- 
ity was sensitive, we still found some resistant individuals, when 
interpreted in the context of Lederberg's experiment, would favor 
the directed (epi)mutation hypothesis. If we had selection of a 
pre-existing mutation or resistance phenotype, then we would 
have expected that all, or no, worms from a same genotype 
would present the same phenotype (resistance or sensitivity). 
From our results, we propose the hypothesis that epigenet- 
ics could be part of the mechanism underlying the resistance 
induction. 

Our ChlP-Seq experiment provided evidences of epigenetic 
variation between resistant and control worms. Our analysis high- 
lighted at least 64 chromatin structure changes between the two 
conditions, with an enrichment profiles in genes implicated in 
catabolic, detoxification, and transfer activity (Two examples can 
be seen in Supplementary Figure 1). No specific pathway linking 
these genes together arose, but it makes sense in the context of a 
response and the subsequent elimination to a drug to see these 
kinds of biological processes being involved (as stated before, 
hycanthone is known to be eliminated from resistant worms a 
certain amount of time after its penetration, while it is stuck 
in sensitive worms metabolism). It is therefore conceivable that 
the drug resistance observed in the surviving worms is due to 
epistatic interaction of multiple gene products rather than an 
(epi)mutation in a single locus. Another interesting finding is 
the presence of a differential H3K27me3 profile in the gene 
body of SmMRPl (Multidrug Resistance Associated Protein 1, 
Smp_171740). In various species, including trematoda, this gene 
is part of a family of transporters playing a role in exclusion 
and elimination of xenobiotic compounds and metabolic toxins 
(Kasinathan et al., 2010). In schistosomes, it was demonstrated 
that SmMRPl has higher levels of expressions in response to 
the widely used schistosomicide praziquantel, and possibly other 
drugs (Kasinathan et al., 2010; Greenberg, 2013). We found a 
significant enrichment in H3K27me3 around exon 21 of this 
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SmMRPl in the control condition. Enrichment of this mark are 
usually associated with increases of facultative heterochromatin 
and gene repression (Trojer and Reinberg, 2007; Bannister and 
Kouzarides, 2011). Our current knowledge of the functional role 
of histone modifications in gene expression in S. mansoni does not 
allow us to affirm that the lower presence of this mark in the resis- 
tant worms would lead to an increased expression of this gene, but 
it stays a promising candidate for further experiments. 

Fourty-seven percent of the genome of S. mansoni is composed 
of repeats (Lepesant et al., 2012). The function of these repeats is 
completely unexplored and, in general, they are not included in 
the analysis of molecular basis of phenotypic variation. We have 
shown earlier (Lepesant et al., 2012), that repeats can be of impor- 
tance in schistosomes, and we decided to study their chromatin 
structure changes upon hycanthone treatment. Only a very small 
fraction of the repeats (0.4%) change their histone modification 
status. In contrast to other organisms, stress (in this case hycan- 
thone treatment) does not seem to lead to an euchromatization 
and potentially mobilization of repeats but has the inverse effect 
(decrease of H3K9ac and increase in H3K27me3). Interestingly, 
the repeats in which H3K27me3 status changes are located in 
satellite-type blocks on specific locations on the genome. It is 
conceivable that they form knobs that separate genome domains 
and influence indirectly gene expression. Further work will be 
needed to understand the impact of these repeated regions on the 
definition of the chromatin structure in which they are included. 

While this ChlP-Seq experiment brings new insight on epi- 
genetic mechanisms and regulation in S. mansoni, it does have 
some limitations. In order to be as faithful as possible to lansma's 
work to reproduce his Type I induction, we exposed the worms 
to hycanthone through an injection in the rodent host. However, 
this generates a degree of variability from one mouse to another, 
as each individual has a different metabolism, immune system, 
clearance rate, and other physiological factor influencing the 
exposure of the worms to the drug. Also, although cercariae issued 
from a same initial miracidium are theoretically clonal, genomic 
instability such as somatic mutation and recombination can occur 
during asexual life stages like sporocysts, as shown by Vieira et al. 
(1991), and Bayne and Grevelding (2003), and their impact on 
the phenotype should be taken into account. Nonetheless, we 
managed to recover worms that survived the hycanthone treat- 
ment at two occasions (initial drug screening and Lederbergs-like 
experiment), from a low genetic diversity strain and without 
any evidences that the gene uncovered by Valentim et al. (2013) 
was involved. This, in addition with the differences of chromatin 
structure observed between surviving and control worms lead 
us to think that induction of resistance in S. mansoni is not 
a myth after all. Nevertheless, we do not think that induction 
should be seen as a specific response of the epigenome to a cer- 
tain environmental condition leading to over/under-expression 
of a specific gene. Instead, this induction phenomenon seems 
to come from the capacity of some strains, or some individu- 
als within a strain, to extend their degree of plasticity beyond 
the genetically defined "default" state through chromatin struc- 
ture changes. This allows them to overcome an environmental 
stress in the absence of a pre-existing stress-response phenotype. 
Depending on the length of the stress, we can imagine that 



the progeny can more or less efficiently inherit the epigenetic 
condition responsible for the phenotype until it is fixed in the 
population (through a DNA mutation for example). This fits well 
with theoretical framework developed recently in which herita- 
ble phenotypic variation is based on at least two components: 
the low-fidelity (epigenetic) system and the high-fidelity (genetic) 
system (Klironomos et al., 2013). In this model, upon changes in 
the environment, the low-fidelity system can generate new phe- 
notypes that can explore the adaptive landscape. The phenotype 
is decoupled from the genotype and appears before the adap- 
tive genotype. The low fidelity of the epigenetic system means 
the adaptative phenotype could theoretically be loss as quickly 
as it appears, unless a strong selective pressure is applied. This 
is consistent with the observations of the metastability of hycan- 
thone resistance in absence of drug pressure (Jansma et al, 1977; 
Brindley et al, 1989), and the stability of the phenotype after 3-5 
successive generations of schistosomes where under drug pres- 
sure (Cioli and Pica-Mattoccia, 1984; Dias and Olivier, 1985). 
Basically, epigenetic variation allows to "buy time" for a few gen- 
eration (in this case, 3-5) until a stable, genetically transmitted 
phenotype occurs. To understand better the underlying mech- 
anism, measuring the expression level of Smp_089320 and of 
candidate genes identified here would be informative in future 
studies. 

Since hycanthone is not anymore used for treatment of bil- 
harzia and the use of oxamniquine is depreciated (Galdino da 
Rocha Pitta et al., 2013), mechanism of resistance is more of 
academic interest. However, experience with hycanthone and 
oxamniquine shows that a new conceptual framework is neces- 
sary to design experiments that distinguish the part of genetics 
and epigenetics in resistance phenotypes. We believe that a stan- 
dard procedure would be beneficial for the testing of resistance to 
current drugs, especially praziquantel, the only molecule actually 
used for mass treatment of schistosomiasis. Resistance to praz- 
iquantel shares striking similarities with early observations on 
hycanthone/oxamniquine. It can be induced through subcurative, 
but increasing doses of praziquantel over six-seven generations, 
and genotyping of two populations of schistosomes showed that 
there was no significant difference in the genotypes before and 
after drug treatment, meaning no selection of specific genotypes 
by the drug (Blanton et al, 201 1; Greenberg, 2013). We think that 
a standard procedure for the next drugs should be similar to the 
one we used to detect resistant individuals from a clonal pop- 
ulation. Mono-miracidial infection would be a good start, as it 
minimizes the genetic variation between individuals and allows 
to detect changes in the epigenotype. A higher number of mono- 
miracidial infections, using several strains of S. mansoni (as the 
induction seems to be strain dependent), and a larger number 
of clones (adult or schistosomula) from each mono-miracidial 
infection should be tested, as we still do not have a clear idea of the 
frequency of the phenomenon. Induction step can be done in vivo 
(for any of the tree types) but screening of resistant progeny has 
to be performed in vitro, as described by Coles and Bruce (1987) 
and Pica-Mattoccia et al. (1992b). This is mandatory to make 
sure that all parasites have the same exposure to the drug. To 
look for heritability, resistant worms can be surgically reinserted 
into new rodent hosts and allow to reproduce. It is critical to 
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have a standardized procedure to maintain the life cycle of the 
parasite (i.e., number of miracidia and cercariae, as well as their 
genotypes, used to infect the hosts) as it is difficult, if not impos- 
sible, in some older publications to find out if variation in the 
numbers of perfused worms come from effect of the drug or 
bias in infections. Using the same intermediate and final host is 
also important as they have an impact on the phenotype and life 
traits of the parasite. In the articles we reviewed, most experi- 
ments were done using the same final host (swiss albino mice), 
but sometimes hamsters were used, too. The geographical ori- 
gin of the mollusk host is not always explicitly stated, and there 
are numerous evidences that it plays an important role in the 
parasite infection success and development (Theron et al., 1997). 
Although this was never explicitly discussed, it could be a cause to 
the contradictory results on hycanthone induction. Finally, both 
genomes and epigenomes of resistant worms would have to be 
compared to the ones of sensitive worms to figure out the ori- 
gin of the phenomenon. Also, from our ChlP-Seq results, we 
saw that it is essential to keep a global, genome scale approach 
including repetitive sequences. The epimutations we found were 
distributed across the genome and are involved in various path- 
ways. A candidate gene approach would probably have missed 
some of the variations we detected. Rather than only focusing on 
find a single gene responsible for the phenotype, we believe that 
the approach we describe would lead to a more complete answer 
on how schistosomes develop resistances. 
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